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Abstract 

We use the nested loop approach to investigate loop models on random planar 
maps where the domains delimited by the loops are given two alternating 



, colors, which can be assigned different local weights, hence allowing for an 

explicit Z2 domain symmetry breaking. Each loop receives a non local weight 
n, as well as a local bending energy which controls loop turns. By a standard 
cluster construction that we review, the Q = Potts model on general random 
CN ■ maps is mapped to a particular instance of this problem with domain-non- 

symmetric weights. We derive in full generality a set of coupled functional 
relations for a pair of generating series which encode the enumeration of loop 
^ ■ configurations on maps with a boundary of a given color, and solve it by 

■ extending well-known complex analytic techniques. In the case where loops 

are fully-packed, we analyze in details the phase diagram of the model and 
derive exact equations for the position of its non-generic critical points. In 
particular, we underline that the critical Potts model on general random maps 
is not self-dual whenever Q ^ 1. In a model with domain-symmetric weights, 
we also show the possibility of a spontaneous domain symmetry breaking driven 
by the bending energy. 



1. Introduction 

The study of random planar maps, which are proper embeddings of graphs in the 
two-dimensional sphere, is a fundamental issue in the fields of combinatorics ^ , as well 
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as theoretical physics . When deahng with the hmit of large maps, many ensembles of 
random maps display the same asymptotic statistical properties, hence define the same 
universality class. Such universality classes are characterized by a number of critical 
exponents, and in some cases, may be entirely determined in terms of a single, properly 
defined, probabilistic object. The simplest such object is the so-called Brownian map 
0], which, in the language of theoretical physics, characterizes the universality class of 
so-called pure 2D quantum gravity and describes in particular the scaling limit of large 
random maps with bounded face degrees. Other universality classes are reached when 
the maps have unbounded and properly tuned face degrees, or when they are equipped 
with additional statistical models, driven to appropriate critical points [Q. 

A particularly important class of statistical models on random maps is formed by 
the 0{n) loop models, which consist in having maps endowed with configurations of self- 
and mutually-avoiding loops. In these models, a non-local weight n is assigned to the 
loops, in addition to a number of local weights and, for < n < 2, the maps are known to 
display non-trivial scaling limits in several regions of parameters 0,^. Phases of dense 
loops and dilute loops have been identified, each characterized by specific non-trivial 
exponents depending on n. 

In 1^, a new approach was introduced to study loop models, which consists in using 
the nested structure of the loops to transform the model into a problem of maps without 
loops, but with faces of large degrees, controlled in a self-consistent way. The advantages 
of the method is twofold: first, by establishing a direct relation between loop models and 
models of maps with controlled face degrees, it allows to identify the universality class 
of loop models in their dense and dilute phases with that of random maps with large 
face degrees in their scaling limit. Such universality classes are entirely determined in 
terms of well-defined non-trivial probabilistic objects, themselves coded by stable trees 
1^. The second advantage of the nested loop approach is that it provides an elegant 
derivation of a number of functional equations for the generating functions of the loop 
model by simply expressing the self-consistency of the mapping to maps with controlled 
face degrees. These equations were already known in some cases from a matrix integral 
formulation of the loop model at hand, and the nested loop approach gives them a clear 
combinatorial interpretation. We applied this method in to address a number of 
0{n) loop models, incorporating in particular some bending energy for the loops. Very 
explicit results were obtained for models where the loops visit only triangular faces, 
which extend the particular results of for loops without bending energy. 

So far, the method was applied only to loop models of the 0(n)-type, which are 
"domain-symmetric" in the sense that there is no parameter in the model that would 
introduce some distinction between the exterior and the interior of the loops. The 
purpose of this paper is to study, on the contrary, loop models on random maps with 
domain symmetry breaking. More precisely, let us define a twofold loop configuration on 
a map as a configuration of self- and mutually- avoiding loops visiting the faces of the 
map, together with a choice of bicoloring, say in red and green, for the domains created 
by the loop configuration. The color of domains is required to change upon crossing a 
loop so that, in practice, the color of all domains is entirely fixed by that of a single 
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one. All the vertices of the map, as well as those edges which are not crossed by loops, 
receive the color of the domain they lie in. The loops are nothing but the domain walls 
for a Z2 variable which lives on the vertices and specifies their color. The notion of 
domain symmetry breaking is equivalent to that of shading for loop models studied in 
||T0[| by random matrix techniques, and in the context of planar algebras. 

In the 0(n)-type models, red and green domains can always be introduced but, 
if so, they play symmetric roles, hence the color is an irrelevant variable. Here we 
shall be interested instead in situations where the red and green colors are treated non- 
symmetrically, either because the parameters of the model are themselves not color- 
symmetric (explicit domain symmetry breaking), or because a spontaneous domain 
symmetry breaking occurs in some a priori color-symmetric model. As we shall explain 
below, the archetype of a loop model with domain symmetry breaking is obtained as the 
loop representation of the Potts model on general random maps (or, equivalently, the 
random-cluster model on general random maps). Besides, more general twofold loop 
models may be considered, including in particular some bending energy for the loops, 
which, as we shall see, may then be responsible for a spontaneous domain symmetry 
breaking in an initially symmetric model. For completeness, let us mention that the 
Potts model on general random maps was previously studied in []TT[| by different methods, 
while its corresponding loop model was also considered in Jiy]. While some of our 
results match those of this latter paper, our approach does not rely on random matrix 
techniques, and we furthermore discuss the critical behaviour of the model. 

The paper is organized as follows: In Section 2, we introduce precisely the twofold 
loop models that we shall study. We show in particular in Section 2.1 how to view the 
Potts model on general random maps as a particular instance of twofold loop model on 
random triangulations. Section 3 describes the nested loop approach to these models via 
a gasket decomposition (Section 3.1), leading to coupled functional equations for their 
generating series (Section 3.2). The solution of these equations is discussed in details in 
Section 4 in a general setting, by introducing an appropriate elliptic parametrization. 
This framework is then used in Section 5 to address the more specific question of non- 
generic critical points, where the maps display asymptotic behaviors which do not fall 
in the universality class of pure gravity. Explicit results are given in the case of the 
Potts model on random maps, as well as for fully-packed loop models with bending 
energy. We show in particular in Section 5.5 how a spontaneous symmetry breaking 
may be generated by the bending energy in some a priori domain-symmetric model. We 
conclude in Section 6. 



2. Twofold loop models 

2.1. The Potts model on random planar maps as a twofold loop model 

As it is well-known, the Potts model on a lattice may be expressed, very generally, 
as a model of loops ||l^-[l^], which turns out to be a twofold loop model as we defined it. 



In this Section, we consider more precisely the Potts model defined on general random 
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planar maps and recall how to rephrase it as a twofold model of fully-packed loops on 
planar triangulations. 

Given an arbitrary planar map M , and for a given integer Q > 1 , we define a Potts 
model on this map by assigning to each vertex v of the map a Potts variable <j{v) taking 
its value in the set {1,2,..., Q}. Given such a Potts variable configuration, each edge 
e of the map is then weighted by t if the two Potts variables at its endpoints v{e) and 
v'{e) have different values, and texp(i^) if they have the same value. In other words, if 
we denote by E, V and F the set of edges, vertices and faces of A^, the weight of the 
Potts configuration {cT{v)}y^v reads 

w{{a{v)}yev) = t'^' n i^^'y{v{e))a{v'{e})), (2.1) 

eeE 

and the partition function Zpotts{-M;t, K) of the Potts model on the map M. is obtained 
by summing over all Potts configurations, namely 

-^Potts vev)- (2.2) 

{cr(-!;)}„gv 

We may eventually consider the grand-partition function 

Zpotts(t,i^;/^.,/^/) = ^^l^^^^f^Zpous{M,t,v), (2.3) 

maps 

where we also sum over all map configurations, with an additional weight fiy per vertex 
and Hf per face of the map. 

Writing exp {Kdaa') = 1 + JSaa', with 

J = exp(K) - 1, (2.4) 

each factor in (|2.1| ) gives rise to two terms (1 or J5o-(u(e))cr(i;'(e))) ^■rid, upon expanding 
the product, we may express the configuration weight as a sum over subsets S of E 
corresponding to those edges for which the term J5cr(„(e))o-(«'(e)) was chosen. This allows 
to recast the partition function of the Potts model on the map Ai as 

Zpotts(yW,i,i^ = ln(l + J))=tl^l 5^ jl^lQ^(^), (2.5) 

SCE 

where the term Q"^^^^ arises from the sum over the Potts variables. Here c{S) denotes 
the number of clusters, i.e. the number of connected components of the graph with 
vertex set V and edge set S (note that a cluster may be reduced to a single vertex). 
In this form, the model may be extended to arbitrary ( no n- necessarily integer) values 
of Q, and is then often called the random- duster model. The random-cluster model on 
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(a) (b) 




Fig. 1: The equivalence between tlie random cluster model on arbitrary 
planar maps and the fully-packed twofold loop model on planar triangula- 
tions, see text. 

random arbitrary planar maps may then be reformulated as a fully-packed twofold loop 
model on random planar triangulations as follows. 

Given a map endowed with a cluster configuration S, as displayed on Fig. l-(a), 
we introduce its dual map A^*, with V* and E* its vertex and edge sets {E* being in 
bijection with E and V* in bijection with F), and we define S* as the subset of edges 
in E* whose dual edge is not in S. In Fig. l-(b), we represent simultaneously the maps 
Ai (red) and Ai* (green), the edges in S and S* being those displayed as solid lines. 
Upon linking each vertex of A^* to all the vertices incident to its dual face in Al, we 
obtain a planar quadrangulation, as displayed in blue on Fig. l-(c), which is made of 
squares whose diagonals consist of an edge in E and its dual edge in E*. By definition 
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of S* , each square has exactly one of its diagonals in SU S*. Thus, adding all the edges 
in U 5"* to the quadrangulation, each square is split into a pair of triangles (forming a 
"diamond"), and we obtain a planar triangulation with vertex set VUV* , see Fig. l-(d). 
By construction, each triangle is incident to exactly one edge in S U S* , and we then 
draw an arch linking the middles of the two other edges. As appears on Fig. l-(e), those 
arches (displayed in purple) form a configuration of closed self- and mutually-avoiding 
loops visiting the faces of the triangulation. The loops visit all the triangles of the 
triangulation, hence we say that the loop configuration is fully-packed. Furthermore, 
the loop configuration is naturally a twofold loop configuration, as we defined it, with 
distinguished red domains containing of all the vertices in V and all the edges in S, and 
distinguished green domains, containing all the vertices in V* and all the edges in S*. 
In Fig. l-(f), we have suppressed the coloring of the edges since it is clearly a redundant 
information with the data of the loop configuration (and of the coloring of one vertex). 

Given conversely a planar triangulation endowed with a twofold fully-packed loop 
configuration, we may easily recover the planar map endowed with a cluster configura- 
tion leading to it by following the above construction backwards. Each triangle being 
visited by a loop has exactly one of its edge sides not visited by a loop. The other side of 
the corresponding edge necessarily belongs to a different triangle and is the unique edge 
side of this triangle not visited by a loop. Concatenating the triangles by pairs along 
their unvisited edge side therefore results into a planar quadrangulation, together with 
a marking of one diagonal for each square. Now the bicoloring of the domains created 
by the twofold loop configuration on the triangulation induces a bicoloring in red and 
green of the vertices of the quadrangulation (which is automatically bipartite). This 
allows to distinguish in each square of the quadrangulation its red and green diagonals 
(colored according to the color of their endpoints), one of them being marked. The map 
Ai is then simply recovered by keeping only the red vertices and all the red (marked 
or not) diagonals, while the subset S is simply the set of marked red diagonals. Note 
that exchanging the red and green colors in the twofold loop configuration leads to a 
different, and dual pre-image for the fully-packed configuration, namely {M.*, S*). 

In conclusion, the above construction provides a one-to-one correspondence between 
planar maps with a distinguished edge subset and triangulations endowed with a twofold 
fully-packed loop configuration. 

Let us now see how to express the weight ( p.5|) in the language of the fully-packed 
loop configuration. Clearly, the factor t'^' jl"^', which corresponds to assigning a weight 
Jt per edge in S and a weight t per edge in i? \ in the original map, is recovered by 
assigning, on the triangulation, a weight Jt per red edge (i.e. an edge strictly inside a 
red domain) of the twofold loop configuration (edge in S), and a weight t per green edge 
of the twofold configuration (edge in S*). Note that each triangle has its unvisited edge 
either red or green and we call it red-facing and green-facing accordingly. The above 
weighting corresponds to assigning a weight \/Jt per red-facing triangle (since an edge 
is shared by two triangles) and y/i per green- facing one. As for the weight Q^^*^) , it may 
be recovered as follows: a cluster (connected component) of S made of e edges and v 
vertices yields a red domain which is surrounded by a number of loops equal to 2 -|- e — v. 



6 



Indeed, it forms a planar map with v vertices, e edges and as many faces as surrounding 
loops so that the above expression is nothing but Euler's formula. Summing over all 
clusters, we deduce that the total number of loops is L = 2c{S) + \S\ — \ V\ since all 
vertices of the map belong to a cluster by construction. The weight Q^^^^ is then 
recovered by assigning a weight 



n 



= \fQ per loop, (2.6) 



a weight l/^/Q per edge of S (or equivalently per red- facing triangle), and a 

weight ^/Q per red vertex in the triangulation. Gathering all the triangle weights, we 
find a weight 

h'^^^ = {Jt/ \/Q)^^'^ per red-facing triangle, 

= t^/"^ per green-facing triangle. 

As for the vertex weights in the triangulation, reintroducing the vertex and face weights 
fj,y and fj,f for A4, we get: 

u^^^ — Uy v/Q per red vertex, 

(2) (2.8) 
= A*/ psr green vertex. 

Prom the above discussion, we have the identity 

ZpotUt,K;iiy,iif)= Yl n#i°°P^(/i(^))#"tH;ngir(/i(2))#^T^^^^^^ 

triangulations twofold fully — packed 
T" loop configs. onT" 

/ ^ / (2) \ 4k S^^en 

X in / J" vertices / J vertices 

'(2.9) 

with K = ln(l + J) and n, h^^\ h^'^\ u^^^ and u^'^^ as above. Note that we may avoid 
introducing vertex weights in the triangulation (i.e. take u^^^ — u^"^^ — 1) by assigning, 
in the first place, a weight //^ = per vertex of the map A4, and no face weight 

(/i/ = 1), i.e. give the map ab initio a weight 

(^/g)-|^lZpotts(>^,^,i^)• (2.10) 

Note also that going from the Potts model on the map M. to that on the dual map M* 
simply amounts to exchanging the red an green colors of the twofold loop configuration. 
We deduce the identity 

(^)-I^IZpotts(A^, t,K = ln(l + J)) = (v^)-l^*IZpotts(A^*, t\ K* = ln(l + J*)), 

(2.11) 

provided wc choose t* and J* such that t* = Jij^/Q and t = J*t* / ^/Q (that is to say 
(/t(2))* = and (/i(i))* = /i(2)), or equivalently 

J J* = Q, Jt^ = J*t*^. (2.12) 
7 



In this setting, self-duality is reached upon choosing J = J* = ^/Q and t = t*. 



In conclusion, the Potts model on general random maps is equivalent to a twofold 
fully-packed loop model on random triangulations, with weights ( ^■6| )-(^). For the 



natural choice /U„ = /Uj = 1, we note that the vertex weights u^^and u^"^^ are not 
color-symmetric, i.e. the model presents an explicit domain symmetry breaking. As we 
shall see, this implies that the triangles weig hts and 

are also, in general, not 

identical at a critical point of the Potts models and their value cannot be obtained by 
a simple argument of self-duality. 

Let us insist on the fact that the correspondence holds between the set of all 
planar maps Ai carrying a Potts model on their vertices, and the set of all planar 
triangulations T carrying a twofold fully-packed loop model on their faces. We have 
not found a natural way to assign weights in T which would allow to control precisely 
the degree of the faces in Ai. This prevents up to now the application of our method 
to the Potts model defined on, say, (the vertices of) random triangulations only. Such 



a Potts model has been studied using matrix integrals ||T5| - |T8|1 and more recently by 
combinatorial techniques [jllf . 



2.2. More general twofold loop models 




Fig. 2: Bending energy: a curvature weight a is assigned to each pair of 
consecutive triangles facing the same domain (i.e. two red-facing or two 
green- facing triangles). 

From the above discussion, the Potts model may be viewed as the archetype of 
twofold loop model with domain symmetry breaking, and a particular attention will 
be paid to this model in the following. Nevertheless, we may want to consider more 
general twofold loop models with domain symmetry breaking. The main generalization 
discussed in this paper consists in introducing, as in ||^, some bending energy for the 
loops as follows: the triangles visited by a given loop form a cyclic sequence of red- or 
green-facing triangles. To each pair of successive triangles along this sequence, we attach 
a curvature weight a if these triangles are the same nature (i.e. two red-facing triangles 
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or two green-facing ones) and a weight 1 otherwise (i.e. one red-facing triangle and the 
other green- facing) , see Fig. 2. The curvature weight has no obvious interpretation in 
the Potts model, however in the random-cluster model ( p.5| ), it corresponds to a corner 
interaction: to each corner of the map, we attach a weight a if its two adjacent edges 
are both in S or both m. E\S. 




Fig. 3: Reformulation of the bending energy for twofold loop models on 
triangulations, as an Ising coupling on tetravalent maps (without loops). 



Alternatively, the curvature weight a may be reformulated as an Ising-like coupling 
for spin variables {^^} as follows. We have seen that concatenating the triangles by pairs 
along their unvisited edge side results into a planar quadrangulation with a marked, 
red or green diagonal in each square. Let us suppress this diagonal and replace it 
by a vertex at the center of the square, with the same color, and consider the dual 
map of the quadrangulation, which is a tetravalent map (see Fig. 3). The color of the 
square centers defines an Ising variable (say q = +1 for red and <j = — 1 for green) 
at the vertices of the dual tetravalent map. Now the loops cross all the edges of the 
quadrangulation. If we now follow the cyclic sequence of triangles visited by a loop, 
we immediately see that having two consecutive triangles of the same (red- or green- 
facing) nature corresponds to crossing an edge on the quadrangulation in such a way 
that the loop remains tangent to the dual edge on the tetravalent map so that this dual 
edge connects two vertices carrying the same Ising variable, see Fig. 4-(a). Conversely, 
having two consecutive triangles of different nature corresponds to crossing an edge 
on the quadrangulation in such a way that the loop also crosses the dual edge on the 
tetravalent map so that this dual edge now connects two vertices carrying different Ising 
variable, see Fig. 4-(b). The curvature weight a is equivalent to an Ising-like coupling 
of the form 1 + [a — l)5<j<j' for all pairs of adjacent vertices on the dual tetravalent map. 
As for the the weights h^^^ and h^'^\ they may be thought of as a magnetic-like coupling 
of the form {h^^^)'^6c,,i + (/i*^^'')^5<j,_i for all vertices of the dual tetravalent map. In 
particular, for n = 1 and u^^^ = u^'^\ the twofold fully-packed loop model on random 
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(a) 



Fig. 4: A pair of consecutive triangles on a loop yields an edge in the 
tetravalent map. In case (a), both triangles are green- facing hence con- 
tribute a curvature weight a. In case (b), the triangles are of different 
nature hence no curvature weight is attached. This corresponds to attach- 
ing a weight a to each monochromatic edge in the tetravalent map. 

triangulations with curvature weight is fully equivalent to an Ising model on random 
tetravalent maps. We shall use this result later in this paper. 

As a last generalization of the model, we may finally release the constraint that 
the loop configuration be fully-packed by allowing for faces which are not visited by 
loops. The loops are still viewed as domain walls separating red and green domains 
and all the vertices of the map, as well as all edges not crossed by a loop and all faces 
not visited by a loop may be colored in red or green accordingly. We shall let these 
unvisited faces have arbitrary degree k >1 and weight them by ^r^^-* accordingly if they 

lie in a red domain, or g^^^ if they lie in a green domain, with (fi'^^'*)fc>i and (fi'^^'*)fc>i 
two sequences of non-negative real numbers. The fully-packed model is then recovered 
by setting g^j^^ = g^j^^ = for all k. We also assign as before a weight u^^^ per red vertex 
and w*-^-* per green one. As for the faces visited by the loops, we still demand that they 
be triangles, among which we again distinguish red-facing and green-facing triangles, 
according to the color of their only edge which is not crossed by a loop. The red- facing 
and green-facing triangles receive respective weights h^^^ and h^'^^ as before. Finally, 
each loop receives a weight n, as well as some bending energy, as defined above, with 
curvature weight a. 

In the language of the Potts model, unvisited faces are known to correspond to 
defects or "dilutions". More precisely, let us consider the so-called dilute Potts model 
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where some vertices are vacant, i.e. carry no Potts variable at all. In consequence, the 



weight (|2.1|) is modified by restricting the product to the set of edges E' whose both 
ends belong to the set V of non- vacant vertices. Performing the cluster expansion as 
before, we find that in ( |2.5|) the sum is restricted to the subsets S of E', and c{S) 
should be understood as the number of connected component of the graph {V',S). 
Going to the loop representation, each vacant vertex of degree k yields a trivial red 
domain reduced to a single red vertex surrounded by a loop of length k. Removing 
this trivial domain and its loop from the triangulation, we obtain a green unvisited face 
with degree k. By adapting the counting arguments leading to (pT^)-(p78|), we find that 
the total number of loops in the resulting map is now L = 2c{S) + \S\ — \V'\, hence 
each loop should still be counted with a weight n = \fQ, while the weights (p.7] )-( p78|) 
are unchanged provided that each green unvisited face of degree k receives a weight 
g^^^ = t'^/'^jly. Here jly is a weight per vacant vertex in the original map, which we can 
control independently from the weight |U„ for non- vacant vertices. As for red unvisited 
faces, they correspond in the Potts language to "dual" defects, namely special faces 
around which all the incident vertices are forced to carry the same Potts variable, see 



for instance [|T9| for a detailed discussion in the context of the regular square lattice. 
Clearly, such a face of degree k yields in the loop representation a trivial green domain 
reduced to single green vertex surrounded by a loop of length k. Removing this trivial 
domain, we obtain a red unvisited face with degree k which should be counted with a 
weight g^j^^ = {Jt)^^'^Q^'^~'^^^^jj.f, where jif is a weight per special face in the original 
map (by convention, edges incident to a special face are weighted Jt instead of (J-M)t). 
We emphasize that, since we are dealing with the Potts model on arbitrary planar maps, 
the degrees of the unvisited faces of both colors are not restricted. While the following 
techniques can be straightforwardly extended to this case to obtain a parametric answer, 
writing down as explicitly as possible the dependence in the initial parameters for the 
dilute Potts model is beyond the scope of this paper. The dilute Potts model on trivalent 
maps was studied in using matrix integral techniques but, as mentioned at the end 
of the previous section, it is unclear how make the connection with our formalism by 
translating the control on the vertex degrees in the loop model. 

In summary, our generalized twofold loop model with curvature weight depends on 
the following parameters: the vertex weights u^^\ u^'^\ the visited triangle weights h^^\ 
h^'^\ the unvisited face weights {g'k^)k>i and ((7fc'^'*)fc>i, the curvature weight a and the 
weight n per loop. 



2.3. Maps with boundaries 

Following the classical map enumeration methodology, we consider the twofold 
loop model on maps with a boundary of length £, i.e. rooted maps (maps with a marked 
oriented edge) whose root face (the face on the left of the root edge) has degree I. 
By convention, we assume that the root face is not visited by loops (note that this is 
always the case for ^ 7^ 3, since the loops only visit triangles in our model). This implies 
that the root face and all its incident edges and vertices belong to a single colored 
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domain, and we say that the boundary is red or green accordingly. For £ > 1, let 
us denote by F^^^ (respectively F^'^^) the generating function for maps with a twofold 
loop configuration and with a red (respectively green) boundary of length i, where for 

convenience we decide to assign to the root face a weight 1 (instead of g^^^ or g[^^)- We 

also set Fq^^ = u^^^ and Fq^^ = u*-^^ for the "vertex-map" reduced to a single red or 
green vertex respectively. Note that, for i = 2, the bivalent root face may be collapsed 
into a single root edge, thus ^2^'* and ^2^'' count twofold loop configurations on rooted 
maps (without boundary), where the superscript indicates the color of the root edge. 

Interpreting the quantities F^^"^ and F^"^"^ as observables in the Potts model involves 
some subtleties. In view of the discussion of the previous subsection, a red boundary 
of length i corresponds in the Potts model to a root face of degree i around which all 
vertices are forced to carry the same Potts variable, which is akin to imposing fixed 
boundary conditions. On the other hand, a green boundary of length i corresponds 
to a vacant root vertex of degree which is close (but not precisely equivalent) to a 
root face of degree i with free boundary conditions. Note however that the quantity 
^2^'* + -F2^'' is equal to the partition function of the Potts model on rooted maps. 

(i) 

In the following, we will obtain a set of equations determining the F^ , which 
will be conveniently rewritten by gathering all these generating functions into the two 
so-called "resolvents" 

£>0 £>0 



3. The nested loop approach to twofold loop models 

We now wish to study the statistical properties of the models introduced in Section 
2, and in particular derive the location and properties of their critical points. To this 
end, we shall use the method introduced in and developed in 0, which consists, via 
some bijective decomposition, in deriving a set of functional equations for the resolvents 
of the model, as defined in ( |2.13| ), and then in solving these equations along particular 
(critical) varieties. 

3.1. The gasket decomposition 

Let us consider a map with a general twofold loop model, as defined in Section 
2.2 above, and with a red boundary of some length i. Following 0, we shall define its 
gasket decomposition as follows: for each loop, we may distinguish its exterior, which is 
the domain separated by this loop which contains the boundary, and its interior, which 
is the complementary domain. In particular, we define the outer (resp. inner) contour 
of the loop as the set of edges that are incident to a triangle visited by the loop, and 
lie in the exterior (resp. interior) domain. Note that all the edges of the outer contour 
of a given loop have the same color while all the edges of the inner contour have the 
opposite color. In other words, each loop has a red and a green contour. Moreover, for 
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a map with a red boundary, the outer contour of any outermost loop is red and its inner 
contour green. The gasket decomposition consists in cutting the map along both the 
outer and inner contours of all the outermost loops. This disconnects the map into: 

• its gasket, which is the map spanned by the edges which were exterior to all loops. 
Clearly, all the vertices of the gasket are red vertices, thus weighted by u^^\ As 
for its faces, they are either regular faces (corresponding to faces that laid outside 
all outermost loops) weighted by g^^^ according to their degree /c, and holes, each 
delimited by the outer contour of some outermost loop. If this outer contour has 
length k (with /c > 1), then the corresponding hole has degree k. 

• for each outermost loop: 

- its ring formed by the triangles visited by the loop. This ring is made of a 
sequence of k red-facing triangles and k' green-facing ones if the outer and 
inner contours of the loop at hand have respective lengths k and k' . 

- its internal map which is the map spanned by the edges which were interior to 
the loop at hand. This internal map is itself a map with a green boundary of 
length k' (the length of the inner contour of the loop at hand) endowed with 
a twofold loop configuration formed by all the loops that laid inside the loop 
at hand. For k' = 0, the internal map reduces to the (green) vertex-map. 

As explained in , all these different components are rooted objects whose root may 
be canonically obtained from that of the original map by some (somewhat irrelevant) 
prescription. In particular, rings are rooted on their outer contour, or equivalently have 
a marked red-facing triangle: we shaU denote by A^i^^^\h^^\ h^^^) {k > l,k' > 0) 
the generating function for rings formed of k red-facing and k' green-facing triangles 
weighted by /i*-^-* and /i*-^^ respectively, and with a marked red-facing triangle. Clearly, 
to enumerate properly all maps with a red boundary, we simply have to assign to each 
hole of degree k a weight 

"E4:?"(''"*.'''")'-i'' (3-1) 

k'>0 

accounting for all possible rings with fixed outer length k and varying inner length k' 
and all internal maps with a green boundary of compatible length k' , together with a 
weight n for the outermost loop that gave rise to the hole at hand. 

3.2. The fixed-point conditions and the functional equations 

We may use the above gasket decomposition to relate the generating function F^^-* 
of maps with a twofold loop configuration and with a red boundary of length i to the 
generating function for planar maps without loops but with a control on the degree of 
their faces. More precisely, let us introduce the generating function J-'^(w; gi,g2^- ■ ■) 
for maps (with no loops) with a boundary of length £ and weight u per vertex and 
gk per face of degree /c, with {gk)k>i a sequence of non-negative real weights (we set 
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J-'o{u; (727 ■ • •) = u hy convention). Then, from the gasket decomposition, we imme- 
diately see that F^^^ is obtained by taking J^£{u; (71, (72, • • ■) at m = u^^^ and at particular 
values of gk, given implicitly by 

9. = S<"+n5;4y^'(h('),h«)Ff. (3.2) 

fc'>0 

Here the two terms come from regular faces of degree k in the gasket and holes of degree 
k respectively. Clearly, we obtain symmetrically by taking J^i{u; gi, §2, ■ ■ ■) now at 
u = u^'^^ and at values of gk given by 

■„ = sf+n 5: feM)F<", (3.3) 

fc'>0 

where 

<,-r^)(/.(^),/.(^)) = <-^)(/.(^),/.(^)). (3.4) 

To summarize, we may write the two identifications 

) = Tdu'^^ ; . . .), = ; G^^, , . . .), (3.5) 

where G^j^^ and G^j^^ are functions of all the parameters of the loop model, implicitly 
determined by the "fixed-point conditions" 

Gl" = si" + n J2 4:-^'(ft<'), fe(^>)^.,(„P);Gf' , Gf , . . .), 

k'>0 

Gf = + n ± ^<y"(ft'", ftM)^.,(uW; G<'>, G<", . . .). 

fc'>0 

It is well-known (see Section 6 in for a proof) that the resolvent 

W(.) = -^'"'fm"--' (3.7) 

A;>0 

may be analytically continued into a function which is holomorphic on the complex 
plane except on a real interval [7-, 7+], with 7_|_ > |7-|, where it has a discontinuity. 
The discontinuity of >V(x) on this cut is given by the so-called "spectral density" 

p(x) = mx--iO)^-W(x + .0)^ j3 

which must be positive on ]7_, 7_|_[ and must vanish at x = 7-1-. Moreover, >V(x) satisfies 
the fundamental relation 

Vxe]7-,7+[, >V(x + iO) + W(a;-iO) = y'(a;), (3.9) 

14 



(3.6) 



where we have set 

V\x)=x-Y,9kX^~^- (3.10) 

fc>i 

Remarkably enough, >V(x), as weU as 7_ and 7_(_, are then entirely fixed by those 
analytic properties and the behavior 

II 

W{x) ~ - (3.11) 

X 

at large x. From the identification (|375|), we deduce that W^^\x) (resp. W^'^\x)) can 
be analytically continued into a function holomorphic on the complex plane except on 
a real interval [7i^'*, 7^^"*] (resp. [7!^'', 74.^'']), with 7^^"* > |7i^''| (resp. 7^'' > |7i^''|). The 
function W^^\x) satisfies, for x G]7i^'', 7^^''[, 

W^^\x + iQ) + W^^\x-iQ) = {V^^^)'{x)-nY^ ^ A^^^,''\h^^\h^^^)F'^^ x^'^ , (3.12) 

k>l k'>0 

and the function W^'^\x) satisfies, for x 

k>l k'>0 

where we have set 



(y«)'(x)=x-5^^7iV-\ z = l,2. (3.14) 

k>i 



The double sum in ( |3.12|) may be rewritten as a contour integral around the cut of W^'^\ 
namely 

k>lk'>0 "^[T- 'T+ ] 

(3.15) 

upon introducing the grand-canonical ring partition function 

A(^-2)(x,y;/.(^),/.(^)) = J2Y1 4]^'\h('\h^'W'- (3.16) 

fc>l fc'>0 

Note that ]7i^'', 7^^ [ must be included in the disk of convergence of y 1— )■ (x, y; h^^\ h^'^^) 

for all X in ]7i^\7^^'*[ in order for the double sum to converge. 

The grand-canonical ring partition function may easily be computed along the same 
lines as in with the result 

xs'i{x) ^ XSi{x) 



y-si{x) ' 2s[{x)' 

(3.17) 
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where we have set 



In the first hne of (|3.17|) , the first term corresponds to k' = in (|3.16|) , in which case the 



ring is made of red-facing triangles only, and the second term corresponds to k' > in 
which case the sequence of red- and green facing triangles may be viewed as an alternance 
of blocks made of red- facing triangles only and blocks made of green-facing triangles only. 
Taking the — In(-) generates cyclic such sequences, and the xd/dx operator amounts to 
rooting the outer contour of the ring, as required. The expression in the second line 
of ( |3.17|) displays both the pole of y i— )■ A^^^'^\x,y; h^^\h^'^^) at y — si{x) and its 



limit xsi{x)/{2s'i{x)) when y — )• oo. In particular, its radius of convergence is 

which implies that si{]j^\'y^\) and must be disjoint in order for ( |3.15| ) to 

converge. Upon exchanging h^^^ and h^'^\ we find a similar expression 

A(^-i)(x, h('\ = + (3.19) 

y-S2{x) 2s^(x) 

with 

^'^"^^ = ahW + il-a^)hWh(^)x ^^-^^^ 

Note that the mappings x i— )■ si{x) and x i— )■ S2{x) are locally decreasing, and they are 
reciprocal 

Si{s2{x)) = S2{si{x)) = X. (3.21) 

The contour integral in ( B.15| ) is easily evaluated by the residue theorem. Its integrand 
^(1^2)^^^ /),(2))^(2)^^) j^a^g a pQig at y = si{x) with residue —xs'i{x)W^'^\si{x)) 

and a residue at infinity equal to •u*^^) xs'i{x)/{2si{x)) since W^'^\y) ~ u^'^^ /y at large 
y. We may eventually rewrite Eqn. ( |3.12| ) as 

T^(^)(x + iO) + Vr(^)(x - iO) -ns[(x)W^''\si(x)) = (V^^h'(x) -nu'^^^^^^, (3.22) 

2si{x) 



and Eqn. (|3.13[) as 



W'^^\x + iO) + W'^^^x - iO) -ns'2(x)W^^\s2(x)) = (V'^^h'(x) - nw^^^^l^. (3.23) 

2s^(x) 

Again W'^'^\x) and W^'^\x), as well as 'j^^ and 7^^'' are entirely fixed by their analytic 
properties, relations ( |3.22| )- p.23D and by the behavior at large x 



,,(1) ,,(2) 

W(^\x)^^, W^^\x)^^. (3.24) 

X X 
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As before, the spectral densities 



2m 

must be positive on ]7i*'',7^*''[ respectively, and vanish respectively at x = 7^''. 

To conclude, let us mention that the functional equations ( p.22|) and ( p.23|) match 



Eqn. (27) of [|T0[ up to notations, in the particular case corresponding to the Potts model 
(a = 1, {V^'^y{x)=x). 



4. General solution of the functional equations 

In this section, we take the functional relations (|3.22| )-( |3r23| ) as starting points, 
and we solve them by adapting the technology of Q , itself inspired from []6|, pmpTl , |2^ . 



The solution consists of four steps: (i) a reduction to homogeneous linear functional 
equations, (ii) a reformulation of these equations using some elliptic parametrization, 
leading to (iii) an explicit solution depending on the position of the discontinuities, 
considered as free parameters. The latter are eventually fixed by matching the divergent 
behaviors so as to fulfill ( |3.24| ). This last step becomes rather technical if one wishes to 
obtain explicit expressions for the position of the discontinuities. We instead focus on 
the construction of the phase diagram of the model by giving (as analytically explicit 
as possible) the location of its critical varieties. This program is achieved in Section 5 
in the case of fully-packed loops. 

We assume that < n < 2, which is the most interesting range for applications. 
This includes the limit n — t- where a single loop survives, and which is different from 
the trivial case 71 = 0. Our analysis holds also for — 2 < n < but we do not discuss it 
here. When |?t,| > 2, only the last step (matching the divergent behaviors) differs. The 
cases n = ±2 are special and could be treated by adapting the remarks of to our 
setting. 

4.1. Reduction to homogeneous linear functional equations 

We may rewrite Eqn. ( |3.22| ) in a more symmetric way upon defining: 

W^^'\x) = s[{x)W^^\si{x)), (V^^'^Yix) = s[{x) (V^^^Yisiix)). (4.1) 
Then we have, for x 

W^^Hx + iO)+W^^Hx-iO)-nW^'''\x) = (y(i))'(x) ^4 2) 



2s; (x) 



and, for x G S2{]'y-\'y+\) 



W^^'\x + iO) + W^^'\x - iO) - nW^^\x) = (y(2-))'(x) + nu^^^^^^ 

2s]^[x) 



(4.3) 
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Here, we used the reciprocity of si and S2, which imphes 

The condition ( p.24|) on W^"^^ translates into the asymptotic behavior 

(2) 

M/(2«)(x)= x^S2(oo). (4.5) 

X - S2(00) 

Since n ^ ±2, we can write 

W^'\x) = Wl2ti^)+W^'\x), W^'^Hx) = W^l:l{x)+W^'\x), (4.6) 
where VFpart(^) ^part(^) particular solutions of ( |4.2[ ), namely 



(1) _ 2(yW)^(x) + n(y(^^))^(x) n(rmW-2^(^)) g;^(x) 
P^'^'^^''" 4-n2 ^ 4-^2 2s; (x) 

(2.). _ n{V^^^y{x) + 2{V(^'^y{x) n(2^(i)-W2)) 
P^'-*^'^''" 4-7i2 ^ 4-n2 2s; (x) 

Eqns. ([4.2|) and ([4. 3)) translate into the relations 



(4.7) 



VT^^^ (x + iO) + W^^'^ (x - iO) - nW^^^ (x) = 0, x G]7i^^ , 7^.^^ [, 
W^^\x + iO) +W^^\x - iO) ~ nW^^\x) = 0, x G S2(]7i^\ 7^.^^[), 



(4.8) 



while Eqns. ( |3.24| ) and (|4.5|) dictate the singularities of the W^^^'s outside their discon- 
tinuity locus. We distinguish two cases: a = 1 and a 7^ 1. 

• a = 1. In this case, si is an affine map, so {V^^^)'{x) and {V^'^^^)'{x) are both poly- 

omials. T. 

for X — )■ 00, 



nomials. Then, the W 's have a (multiple) pole at 00, and since s'/ = we find that. 



(4.9) 



T77(i). X 2(y(i))'(x) + n(y(2^))'(x) /I 

VT^ ^(x) = ^ ^ + + o - 

4 — X \x 

_(2) ^ n(yW)-(x) + 2(y(2s))/(^) ^ ^ /I 

4 — X \x 

a 7^ 1. The pole of si is now located at the finite value 

«2(oo) = ^^-^1^ ^ 00, (4.10) 
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so (V^^''y(x) is a polynomial and (V^'^^^y(x) is a rational function with a pole at 82(00) 

(and no r( 
and since 



(and no residue at a; = 00). Then, the W 's have (multiple) poles at 00 and 82(00), 



s'lix) -1 



2s[{x) X — 82(00) ' 

we find that, when a: — >■ 00, 

^„ ^ _2(^<")'(x) ^ 2(2ua)-n„(^)) 1^/1 

a: \x 

—(2),, n(V(i))'(a;) n(2«(i) - nw^^)) 1 /i 

W '{x) = \ ^ + ^ — 2 ^ - + o - 

4 - 77,2 4 -n2 a; 

while, for x — >■ 82(00) 



(4.11) 



(4.12) 



-d) _ n(V(^-))-(x) n(n^(^) - 2^(^)) 1 ^ 
W- 4_^2 + 4_^2 a; - 82(00) ^ 

^(2) ^ _2(y(^ ^ 2^^)-!^ 1^ ^ 
4-n2 4-n2 a; - 82(00) 



(4.13) 



4-2. Discontinuities 



We shall consider for the time being that the discontinuity edges 7^'' are given. 
They will have to be determined later in terms of the parameters of the model (see 
Section 4.6). As explained in Section 3.2, ]7i^\7^^^[ and S2(]7i^\ 7+ [) do not overlap. 
In particular, when the parameters are small enough, [7!^'', 7^^'*] and [7^^^ 7^"*] are small 
neighborhoods of and, since 82(0) = l/{ah^^^) > 0, 82 ([7!^'', 7+^]) is a small interval 

away from 0, and to the right of [7i^\7^^'']. Since 82 is locally decreasing, with a pole 
at 

we have more precisely, for small enough parameters, 

S2(7i'^) > ^2(7^) > 7i'^ > > ^^y. (4.15) 

(2) (2) (2) 
When the parameters varies, S2([7_ ,7_|_ ]) may swallow 00 and 82(7- ) jumps to the 

interval ] — oo,7^^'*[. It should become clear later that the solution is perfectly regular 
across this jump. 

In any case, the real line can be decomposed in four sets with non-overlapping 
interiors: 

R U {00} = C(^) U C(2) U D^^^ U D^'^\ (4.16) 

where: 

CW = [7i^7i^^ C(^) = 82([7i^\7f]), i^« = [7i^\«2(7?^)], (4.17) 

and D^^^ is swept by following RU {00} in positive direction starting from 82(7!^^) and 
^ (1) 

arriving at 71 . 
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Fig. 5: Mapping of the complex plane with two cuts onto a rectangle via 
Eqn. (4.18). 

4.3. Elliptic parametrization 

Our strategy is very similar to the one used in the context of 0{n) loop models 
with bending energy P], which correspond here to the special case 
is natural to introduce the elliptic integral which parametrizes the complex plane with 
two cuts C^^"^ and C^'^\ namely: 

v{x) = c r (4.18) 

v/±(? - I'-'m - li'^m - s^h'^'mc - .2(7f )) 

We choose the sign ± so that the polynomial under the squareroot is non-negative on 
D^'^\ and a positive constant c enforcing v{D^'^^) = [0, 1], see Fig. 5. Now, if we follow a 

path starting from 52(7!^^) along C*-^^ +iO, v runs over the positive imaginary axis, and 
we introduce T g]0, +00] such that f(C'^^-* + iO) = i[0,T]. Actually, v maps the upper 
half-plane to ]0, l[+i]0, T[, and the lower half-plane to ]0, l[+i] — T, 0[. Finally, we may 
check that v maps 7^^'* to 1 =F i^"- We introduce 

wi =v{+oo + iO), W2 = v{s2{+oo + iO)). (4.19) 

If the local weights in our model are nonnegative, wi and (1 — W2) must belong to 
]0, l[Ui]0, T[. The inverse function x{v) is a priori defined in the rectangle 71 =]0, l[+i] — 
T, T[, called the "physical sheet". By Schwarz reflection principle, it can be analytically 
continued to the whole complex plane as an even, doubly periodic function of periods 2 
and 2iT. Then, 

± (4.20) 

V(^)(.) = = x'(.)(W^^))'(x(.)), 

dv 
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are defined as odd, doubly periodic meromorphic functions witli periods 2 and 2iT. Let 
us define similarly, for v in the physical sheet 71, 

u^'\v) = x'{v)W^"\x{v)). (4.21) 

We are going to extend these functions to the whole complex plane as follows. 

• Since {W^^\x))i=i^2 are continuous across D^^\ we can extend a;*-*-* as a meromor- 
phic function on the strip 5 =]0, l[-t-iR by the relation 

uj^'\v + 2iT) =uj'^'\v). (4.22) 

• Since W^^\x) is continuous across C-^-*, we can extend u^^"* as a meromorphic 
function on the strip S^^^ =] — 1, l[-|-iR by the relation 

uj(^)(-v) = -uj(^)(v), (4.23) 

• Similarly, W (x) is continuous across so we can extend o;*^^-' as a meromorphic 
function on the strip 5'-^-* =]0, 2[+iR by the relation 

uj^^\2-v) = -J^\v). (4.24) 

• The behavior of W^^"^ on its discontinuity (first line of Eqn. ( [4.8|) ) translates into 

V?;Gl + iR u^^\v) + u^^\v - 2) - nu^^\v) = 0, (4.25) 

and this relation allows to extend the definition of cj*^^^ to the strip ]l,2[+iR and 
via ([4.23P to the strip ] — 2, — l[-|-iR so that u^^^ is now defined as a meromorphic 
function on the larger strip (S^^^)' =] — 2,2[-|-iR. 

• The behavior of W on its discontinuity (second line of Eqn. (|4.8|) ) translates into 



V^;eiR uj^^\v) +uj^^\v + 2) -nJ^\v) = 0, (4.26) 

and this relation allows to extend the definition of u^"^^ to the strip ] — 1,0 [+iR, 
and via ( [4.24]) to the strip ]2, 3[+iR so that u^'^^ is now defined as a meromorphic 
function on the larger strip {S^"^^)' =] — 1, 3[+iR. 

By recursion, we can then analytically continue w*^*^ as meromorphic functions on 
the whole complex plane which satisfy, for all f G C, 

uj^^\v + 2iT) =J^\v), 

uji^)(-v) = -J^)(v)^ (4.27) 
uj^^\v) + uj^^\v - 2) - nuj^^\v) = 0, 
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and 

u^^\2-v) = -J^\v), (4.28) 
J^\v) + J^\v + 2) - nuj^^\v) = 0. 

The spectral densities introduced in Eqn. ( |3.25| ) can be computed with the new 
parametrization 

P [Hv + ^)) 2i7rx'(t; + l) , 

(4 29) 

p^^' (si(x(v))) = ,,,,,,,, \ -1^, . 



Building the solution of Eqns. ( U-^lj ) and (U-SSj) 



Let us denote by T the operator of translation by 1 acting on the space of mero- 
morphic, 2iT-periodic functions. The last lines of (^4.27 ) and ( [4.28|) can be combined 
into 

(T2+nT + l)(T2-nT + l)cjW = 0, i = l,2. (4.30) 

Notice that, in the correspondence between the twofold loop model and the Potts model, 
we have n = a/Q, and the equation above does not depend on the sign chosen for the 
square root. We parametrize the weight per loop as 

?i = 2cos7r6, (4.31) 

and the assumption < n < 2 is achieved by taking < 6 < 1/2. Then, we can 
decompose any solution of ( |4.30| ) into a linear combination of pseudo-periodic functions, 
since the corresponding linear operator has the kernel decomposition 

Ker(T2 + nT + 1)(T^ - riT + 1) 

(4 32) 

= Ker(T - e'^^) © Ker(T - e"'^'') © Ker(T + e'""^) © Ker(T + e"'^^). 

There exists a unique function v t— )■ Cb{v) in Ker(T — e^'^^) with a simple pole at f = 
of residue 1 and no other singularity modulo Z + iTZ, i.e. satisfying 

Cb{v + l)=e'-\tiv), Ct{v + 2iT)=aiv), Ctiv) (4.33) 

It can be expressed as a ratio of Jacobi theta function as 

_ Mv + ^bT/2\iT) ^;(0|iT) 

^'^""^ - ' — — ^i(i6T/2|T)- ^^-^^^ 
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This function differs from that used in ||^ by a modular transformation iT — t- — l/(iT). 
Besides, LiouviUe theorem imphes that the only bounded function in Ker(T — e^'^^) is 
zero, so any function in this subspace is uniquely characterized by the behavior at its 
poles in a fundamental domain, and can thus be expressed as a linear combination of 
derivatives — w), where w denotes a pole in a fundamental domain, given that 

such a term yields a contribution 

dla(v-w) = \i+i +0{l) (4.35) 

to the divergent behavior. The description of the other subspaces in (|4.32|) is obtained by 



changing (v, h) to (— f , 6), (f , h) and (— f , h) respectively, where 6 = 1 — 6 is characterized 
by — n = 2 cos nh. 

Now, we have to impose the second line of ( [4.27| ), namely that o;*-^-' be odd. The 
general odd solution of ( ^.3Up with a simple pole at w for some point w in the physical 
sheet TZ takes the form 



Z^^\v] w. A) = A- [Cb{v -w)- Cb{-v - w)] + A+ kb{v + w)- Cb{-v + w)^ 

^ (4.36) 
+ A_ [C-,{v -w)- Cii-v - w)] + A+ [C-,{v + w)- Cii-v + w)] 

with a vector of complex coefficients A = {A-, Aj^, A-, A^). Note that this solution 
also has a simple pole at 1 — w, which also lies in the physical sheet. We may then 
generate the general odd function with an (Z + 1)*^ order pole by differentiating I times 
Z^^\v;w, A) with respect to w. As we shall see in the next subsection, u^^\v) will be 
obtained by taking linear combination of such functions (with several values of w, A 
and I). Note that, according to the third line of ( [4.27| ), a term Z^^\v;w, A) in u}^-^\v) 
gives rise in ui^'^\v) to a term 

Z^'^hv; w,A) = - IZ'^^^v; w, A) + Z^^\v - 2; w, A)l 
n 

= A_ [e-'^%{v -w)- e'^%{-v - w)] + A+ [e-'^'Cb{v + w) - e'^'Cb{~v + 

+ A_{e'^%{v -w)- e-'^%i-v - w)] + i+ [e'^%iv + w) ~ e-'-\~,{-v + w) ] 

(4.37) 

and those terms would precisely be canceled in (^4.28|) . Moreover, those terms give rise 
in p'^^\v) x'{v + 1) to the term 

A^'^\v; w, A) = -^{Z'^^^v - 1; w, A) - Z^^\v + 1; w, A)| 
2i7r 

- sm7r6 |^^ ^^^^^ -w)+ Cbi-v - w)] + A+ [Cbiv + w) + Cbi~v + w)] 

+ i_ [C-,{v -w) + Cd-v - w)] + A+ [C-,{v + w) + C-,{-v + w)] } 

(4.38) 



TT 
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and in p^-^^v) s'i{x{v))x'{v) to the term 

ZITT 

= ^^^^\^A-[Cb{v - w) + Cbi-v - w)] + A+[Cbiv + w) + Cb{-v + w)] 

- i_ [C~,{v -w) + Cd-v - w)] - A+ [C-,{v + w) + C-,{-v + w)]y 

(4.39) 

Notice that A^'^\v; w, A) differ from A^^\v; ty, A) by a simple sign change, namely 

A^'\v; w, A+, i-, i+)) = A(2) (v; w, A+, -A+)). (4.40) 

4.5. Matching the divergent behavior 

We now describe a choice of basis functions which is appropriate for the type of 
divergent behavior encountered in our problem. For convenience, we distinguish the 
case a = 1 from the case a 7^ 1. 

• a ^ 1. For i = 1, 2, let us write the Laurent series of cj*^*-* around Wi as 

u^^){v) = V , +Q(1), (4.41) 

l>0 ^ 

with w\ and W2 as in ( [4.19| ). The coefficients oii\i can be expressed from the first line of 
( [4.12| ) (for z = 1) and the second line of ( |4.13| ) (for i = 2). In particular, we have 

2(2w(i) - nw(2)) 2(nw(i) - 2u^^^) 
= 437^2 ' «2|o = ^3^^2 , (4.42) 

and we note that in the fully-packed case ((yW)'(2:) = x), a.n vanishes for I > 2. 
Furthermore, by the second line of ( |4.12 ) and the first line of ( [1.13| ), we have 



^(3-)(^) = ^^1' +o{l). (4.43) 



We must ensure that w*^*) is regular oX v = \ — w\ and v — \ — W2 iox both i = 1,2. 
This leads us to define vectors Aj G for j = 1, 2 as the unique solutions of the 4x4 
systems of equations 

Z^^\v■w,A^) Z^^\v■,w,A^) = 0(1), 

nl2 ^^-^^^ 
Z^''\v-w,A^) Z^^\v-w,A^) = 0(1), 

v^w V — W v^l — w 



24 



UJ 



where we introduce the differential operators 



The corresponding spectral densities read 



x'{v + 1) 



^ ^ ' s[{x{v))x'{v) 



(4.45) 



and 

v—>-w V — W v^l — w 

We find 

Ai = ^ • (1, 1,1,1), ^2 = \- {e'^\ e-'^\ e-'^\ e'^^). (4.46) 
Now, if we define the two sets of functions (j = 1, 2) 

Zf\v) = Z^'\v;w„A,), Zf\v) = Z^'\v;w„A,), 
= A(i)(^;^^,A,), Af\v)=A'^^Hv;Wj,A,). 

we deduce that the desired solutions (matching all required behaviors) are 



(4.47) 



(4.48) 



(4.50) 



• a = 1. Then, wi = W2- The basis of functions we described for a 7^ 1 is still well- 
defined and appropriate. Nevertheless, since a = 1 has an interest of its own, we prefer to 
define another basis of functions (related to the former one by a linear transformation), 
which takes Eqn. (|4.9| ) as a starting point: let us now write the Laurent series of u^^^ 
around Wi = W2 as 



J^Hv) = J2 f ^^"v+i +^W- (4.51) 



z>o 

In particiilar, we have 



Pilo = -u^'\ P2\o = -u^^^- (4.52) 
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It leads us to define vectors Bj G for j = 1, 2 as tfie unique solutions of the 4x4 
systems 

Z^^\v;w,B,) Z^^\v;w,B^) = 0(1), 



Z(^\v;w,B^) = 0(1), Z(2)(^;;u;,Bi) = 0(1) 



(4.53) 



and 



Z^^\v;w,B2) = 0(1), Z^^\v;w,B2) = 0(1), 

Z(2)(t;;«;,B2) Z^^\v;W, B2) = 0(1). 

v^w V — W v^l—w 



(4.54) 



We find 



Bi = — ^ ■ (e^-^ -e-^-^ -e-^-^ e^-^), B2 = — ^ ■ (-1, 1, 1, -1). (4.55) 

4l SlUTTO 4l SlUTTO 

If we now replace the definitions ( [4.47| ) by 

zW(t;) = zW(t;;«;i,B,), zf{v) = Z'^^\v■,w^,B,), 

(4.56) 

aW(^) = AW(t;;ti;i,B,), Af ^(t;) = A^^H^; B,), 



we deduce that the expressions ( [4.48|) and ( [4.50|) for u^'^^ and p*^*^ are still valid provided 
we redefine the differential operators as 

4.6. Consistency conditions and position of the discontinuities 

We now come to the crucial question of the determination of the cuts in our problem. 
The following discussion is valid both in the case a 7^ 1 with the expressions ( [4.47] )- 
( [4.50|) , and in the case a = 1, with the modified definitions (|4.56| )-( ^757| ). The spectral 
density p^^\x{v + 1)) must vanish at x{v + 1) = (corresponding to = 0) and 

at x{v + 1) = 7^^'' (corresponding to f = — iT). Similarly, p^'^\si{x{v))) must vanish 

(2) (2) 
at x{v) = S2(7_ ) (corresponding to f = 0) and at x{v) = S2{'j\- )) (corresponding 

to f = — iT). This gives four independent conditions which determine in principle the 
positions of 7^!^ By (lOOD , we can decompose A^*'' for j = 1, 2 as 



(2) ~ (^-^S) 
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where Aj only involves the function (i, and Aj only involves C^, see again ( [4 .381) and 
( [4 .SOI) . We deduce decoupled conditions 



PiAi(0)+I?2A2(0) = 0, 

(4.59) 

I)iAi(-iT)+p2A2(-iT) = 0, 

and 

PiAi(0)+I?2A2(0) = 0, 

(4.60) 

I)iAi(-iT)+p2A2(-iT) = 0. 

Last but not least, we must ensure that the spectral densities are positive in the 
interior of their support, which in general selects a unique solution to the above equa- 
tions. Checking global positivity is a difficult task, that we do not carry otherwise than 
graphically in the examples of this article. However, one can often check analytically if 
the densities are locally positive at the edges of their support. We believe that this is 
actually sufficient: 

Conjecture. Spectral densities for enumeration of maps in loop models are locally 
positive at the edges of their support iff they are globally positive on the whole support. 

This conjecture has been verified in all examples (maps with bounded face degree, maps 
with loops) we know of, and is also verified in the examples given later in this article. 



5. Non-generic critical points 

We shall now focus our discussion on the description of non-generic critical points 
(see for a discussion on this terminology), corresponding to values of the parameters 
for which the generating functions for maps with a large boundary behave as 

with an exponent u differing from the non-critical value u = 1/2 and the generic critical 
value ly = 3/2. By transfer, the resolvents W^^\x) should display non-trivial singulari- 
ties of the form (7^"* — xY . This can only happen when 

7i'^ = «2(7f ), (5.2) 

i.e. when T = 00 in the elliptic parametrization. Our main goal will be to derive the 
location of the non-generic critical variety in the parameters of the model. We carry this 
analysis as explicitly as possible in the case of fully-packed models, i.e. gl' ^ and thus 
{V^^^y{x) = {V^'^^y{x) = X. The parameters of the model are then u^^\u^'^\ h^^\ h^'^^ 
the bending energy a and the weight per loop n. We shall put emphasis on: 

• the symmetric model: u^^^ = w*-^-* (and consequently h^^^ = h^"^^). 

• the Potts model in absence of additional vertex or face weights: ^„ = /U/ = 1, so 
that u^^^ = y/Q and u^"^^ = 1 according to Eqn. (|2.8| ). 

• the limit n — t- for general values of the parameters. 



27 



S2(oo) 







X 














i? 


My':') 




















-m 


•- 


• 


> 


h 


9 


m- — 1 


> 1 
















W2 




■ 



Fig. 6: The elliptic parametrization in the limit 7^^'' — ?■ 52(7^'') via 
Eqn. (5.4). 

5.1. Characterization 

When Eqn. ( |5.2|) is satisfied, the elementary function Cb{v) reduces to 

Cb{v) = 2m— — . (5.3) 

It is indeed easy to check that the right-hand side matches the required properties (|4.33| ) . 
Besides, the elliptic parametrization reduces to a trigonometric one, see Fig. 6, 

x(v) = X-6 , (5.4) 

COS TVV — COS TTWi 

where the three parameters wi, X and S are related to the position of the cuts via 

7^^ = ^2(7f) = A-5, 7i^^=A---^^, .2(7i^^) = A. (5.5) 

1 + cos TTWi 

The point x = 7^^^ is now mapped to f = — iT = — ioo, and we have 

x{v) = 7^^^ -25(cos7r«;i - l)e-^''^ (5.6) 

v—^ — ioo 

The value of W2 also fixed in terms of (wi. A, d) since, from Eqn. ( [4.14|) , we have x{w2) — 
S2(cxd), namely 

cos 7rW2 — cos TTZWi (a^ — l)/i'.^J 

We start with three useful observations: 

• wi ranges over ]0, l[UiR;^, and the condition that 7^^"* > 7i^'' implies that 5{1 — 
cosTTWi) remains positive. 
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We note that the asymptotic behaviors 

x'{v) ~ 2i7r 5(cos TTtWi — 1) e' 



x'{v + 1) ~. -2i7r 5(cos ttwi - 1) e-^''^ 



-lOO 



have opposite signs, and that si is locaUy decreasing on the real hne. Hence, the 
relative sign of p^^\x{v + 1)) and p^'^\si{x{v))) near v — ioo is also the relative 
sign of p^^\x{v + l))x'{v + 1) and p^'^\si{x{v))) s'i{x{v)) x'{v). 
• For bookkeeping, we note the asymptotic expansions deduced from (|5.3D 

Cbiv) = 2i7ry e-^"(2-6+2fc).^ ^ -2i7ry e-^"(''+2'=)^ (5.9) 

v—^ — ioo ' ^ >■— icxD < ^ 

/c>0 A;>0 

Among the four conditions ( [4.59| ) and ( [4.60D , the two conditions that p^^^{x{y ^V)) 
and p*^^)(si(x(f ))) vanish at f = —ioo (second lines with T — )■ oo) are better rephrased by 
saying that the coefficient in front of the terms of order e~^^^" (arising from Cb(— f izwj)) 
and Q-^'^i^-^)'" (arising from (^^(— t^i^j)) must vanish, since they would lead to divergent 
terms in the densities (see for a similar argument). This being enforced, the dominant 
order in the densities is q-'^'^'O-^^)'" and arises from terms involving C^^iv ± Wj), but 
Eqn. (|4.40| ) and the observations we just made imply that this term comes with an 



opposite sign in p^^^ and p^'^^ . Thus, the positivity of the densities may only be achieved 
if we impose an extra relation, namely that the coefficient in front of this term vanish. 
Then, the dominant order in the spectral densities when v — )■ —ioo is and 
comes from terms involving ±Wj), exactly as for a non-critical 0(n) model (with 
n = 2 cos irb) in its dense phase. The exponents of our non-generic critical twofold loop 
model are therefore generically the same as those of a usual loop model in its dense 
phase. Extra cancellations may drive these exponents to other values: if the coefficient 
in front of the dominant order Q-^'^i'^-^)'" for one of the spectral densities happens to 
vanish at some values of the parameters of the model, the dominant order for this 
spectral density becomes in practice e~^^(^+^)" (arising from the sub-dominant order in 
Cbi^v ± Wj)), leading to exponents whose values are now characteristic of the dilute 
phase of the 0{n) model. 

To summarize, assuming a non-generic critical behavior imposes 7^^^ = ^2(7^'*), 
together with the "extra relation" that we just discussed. The non-generic critical 
variety is therefore of co dimension 2. 

5.2. Strategy for the fully-packed model 

From now on, we restrict ourselves to the fully-packed model. If the weight per loop 
n and the bending energy a are considered as fixed, we may parametrize the non-generic 
critical variety by expressing h^^^ and h^'^^ on this variety as functions of u^^^ and u^'^\ 
or equivalently as functions of the reduced variables 

u= Vu(^^u(^\ r = ^. (5.10) 



29 



Carrying out this program is rather technical, so let us explain first our strategy. It is 
easily seen from Eqns. ([4.12| )-( ^?T3|) that, in the fully-packed case, uj^'^\v) has a triple 
pole at wi and W2- On top of the residues known by Eqn. ( |4.42| ) (or Eqn. ( [4.52| ) for 
a = 1), there are four unknown Laurent coefficients (Q!i|j)i,j=i,2 (or (/3i|j)i,j=i,2)- Their 
expression in terms of wi, W2 and the parameters of the model can be derived directly 
from Eqn. ( [4.13|) (or Eqn. ( |4.y| )). We prefer to invert these relations to write h^^\ h^"^^ 
and the position of the cuts in terms of wi , W2 and the Laurent coefficients. On the other 
hand, the 4 equations ensuring that the densities vanish at the edges, can be viewed as a 
4x4 linear system determining the 4 Laurent coefficients in terms of wi,W2, u, r. Thus, 
we get expressions for h^^\h^'^^ and the position of the cuts in terms of wi,W2-,u,r. 
Furthermore, using the solution of the 4x4 system to enforce the extra relation^ we 
will obtain a first relation between wi,W2 and r. Enforcing Eqn. ( |5.7| ) gives a second 
relation between wi, W2 and r (for a = 1, it is replaced by the even simpler requirement 
that wi — W2). Hence, wi and W2 are entirely fixed by the value of r. Eventually, we 
check for the local positivity of the densities at 7j_ , which may rule out certain values of 
parameters for the solution, and we check numerically for some admissible values that 
the densities are indeed globally positive. 

In this way, we can arrive at a full analytical solution of the non-generic critical 
variety in the case a = 1, and a semi-analytical one in the case a 1. Although a = 1 
can be obtained by taking the limit a — ?• 1 (which is regular), we prefer to treat this 
case separately since it is less cumbersome than the general case. 

5.3. Fully-packed model without bending energy [a = 1) 

Recall first that a = 1 corresponds to wi = W2. We use the set of functions 
introduced in Eqn. (|4.56| ) and the differential operators of Eqn. (|4.57| ), and solve the 
4x4 system determining the position of the cuts for the Laurent coefficients (/Si|j)i,j=i,2- 
It can be written as 



Vi [sin7r6(l — Wi)~\ + V2 \_s\.mxhwi\ = 0, 
r»i[sin7r(l - 6)(1 - zwi)] - 1^2 [sin7r(l - 6)m;i] = 0, 



cos7r(l - b){l - wi)l 

: "T 

SmTTWi J 

"COS 7r6(l — wi] 



cos 7r(l — b)w 



smnwi 
cos Tcbwi 



smnwi 



(5.11) 



smnwi 

Then, the extra relation reads 

I^i [sin7r(l + b){l- Wi)] - I?2 [sin7r(l + b)wi] = 0, (5.12) 
and if we plug in the solution of Eqn. (|5.11| ), we arrive to the surprisingly simple relation 



Kb{l -Wl)' 



KbiW 



b cosTvbw siuTTw — sinirbw costtw. 



(5.13) 
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which fixes wi as a function of r only. It is easy to see that, when r assumes positive 
values and 6 g]0, l/2[, Eqn. ( |5.13| ) has a unique solution wl{r) g]0, 1[, and no solution 
in iR_(_. Then, we find that the spectral densities behaves near 7_j_ as 



1 /4'u sin7r6\i/2 / sinnwi 



ioo 27T \ b gb{wi) J Vsin7r6(l — wi 



1/2 



1 /4w sin7r6\ 1/2 / sinTTzwi \ 1/2 



v^-ioo 27T \ b gb{wi) / Vsinyrfeiyi 



r{l-b)v 



(5.14) 



and near 7!*^ as 



b Mu sin7r6\ ^/^ 



(:^_^W)l/2 



x^-y'^y V^T^^^9b{wi) ^ (sin7ri(;i)i/4(sin7r6(l - ti;i))3/4' 
b /4w sin7r6\i/4 (x - 7!^^/^ 

,21 hn,Jii)^ \ I ('sin TT-j 



(2) bgb{wi) J {simrwiy/'^{sm7Tbwi)^/'^ 



(5.15) 



Hence the spectral densities are locally positive at the edges of the cuts for all values of 
wi g]0, 1[, and we checked numerically for several values of b and wi that they remain 
globally positive along the cuts. 



a = 1 



0.22 



0.18 




b = OA 



Fig. 7: Non-generic critical values of /i*-^-* and /i*-^-' as functions of r 
for b = 0.4, a = 1, u = 1. 



Next, as explained in Section 5.2, using Eqns. ( ^y9|) - (|4.21|) - (|4.52| ) and the solution 
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of Eqns. (|5.11| )- (|5.12| ), we arrive after some algebra at 



7^^ 



4w sin7r6\ -1/2 1 A/sin7r6(l — wi) smnbwi 



hgh{wi) J h simih ^Jsi^n.'KWl 
4m sin7r6\ -1/2 1 \/sin7r6wi sin7r6(l — wi) 



bgb{wi) J b simrb ^/sm^lWl 
4u sin7r6\ 1/2 b cos7r6(l — Wl)^/sm^lWl 
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bgb{wi) ) A/sin 7r6(l - wx) 

(1) ^Au sin7r6\ 1/2 b sinTriyi cos7r6(l — w\) — (1 — costtzwi) sin7r6(l — 



(5.16) 



(2) 

7+ 



bgbiwi) J y^sinnwi sm7r6(l — Wi] 

4u sin7r6\i/2 b cos nbwiy/ sin tcwi 



bgb{wi) / y/siii Ttbwi 

(2) /4m sin7r6\ 1/2 b cos7r6?x;i sinTrwi — (1 + cosTriOi) sin7r6wi 



bgbiwi) ) a/ sin TTtwi sin -Kbwi 

where we have set 

gb{w) = ^/Kb{wi)nb{i - wi). (5.17) 

This is a parametrization of the non- generic critical variety in terms of u, and wi which 
is itself uniquely determined by r (Eqn. ( [5.13| )). Notice that the ratio h^^^ /h^'^^ has a 
simple expression in terms of wi, namely 



^ / sin7r6wi 
/i(2) y sin7r6(l — zwi) 

For illustration, we have plotted in Fig. 7 the values of h^^^ and h^-^\ as given by (|5.13| ), 
( ^.17|) and ( p.l6|) , at 6 = 0.4, u = 1 and for varying r, i.e. u^^^ = y/r = l/u^'^\ 

• In the limit case n — ?■ (i.e. 6 = 1/2), we obtain the simple relation 

r = tan^{7iwi/2), (5.19) 
so that we may write, in terms of the parameters u, r, 

and the cuts are symmetric with respect to 0: 

7^^ = ±2 u'/^ r^l\ = ±2 v}l^ r-^l\ (5.21) 

• The symmetric case (i.e. r = 1) is particularly simple as Eqn. ( 5.13| ) yields twi = 1/2 
so that, from ( [5.16| ) or ( |5.1ij| ), 

/,(i) = /,(2) = ^ (5.22) 

2V2^v^2Tn 
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Fig. 8: Critical couplings as a function of Q for the Potts model (plain 
line), and for the symmetric model (which correspond to self-dual values, 
dashed line). Here, w = 1 is assumed. 



as already known from |2^. The model is then fully red-green-symmetric. 

• Let us now return to the Potts model with ii^ = fJ-f = 1, i.e. r = 2cos7r6 and 
u = cos7r6. We get from (|5.13D and ( ^.16|) expressions for the critical couplings t 



and J via 

t = (/i(2))2, J= V^(/i(i)//i(2))2, (5.23) 

as obtained by inverting Eqn. ( p.7| ). We have plotted these values as functions of 
Q = 4cos^7r6 in Fig. 8. We note that, for Q 7^ 1, those are not the self-dual values, 
which correspond instead to a symmetric model with u^^^ = u^"^^ = 1, hence h'^^^ = h^'^\ 
These self-dual values for J and t are plotted for comparison in Fig. 8 in dashed lines. 
We insist here on the fact that the Potts model on general random maps is not self-dual 
for Q 7^ 1 (and the value of its critical couplings cannot be obtained by a simple duality 
argument), although its critical point clearly lies in the same universality class as the 
domain-symmetric model. This is in contrast with the Potts model on regular lattices, 
where the self-duality argument is commonly used to determine critical values of the 
couplings. 



5.4- Fully-packed model with bending energy (a 7^ 1) 



7^ 1, we use the Laurent coefficients (ai|j)i,j=i,2- By comparing 
(|4.13 ) and (|4.21 ), we first obtain explicit expressions for h^'^\ 



When a 
Eqns. iWm . 

A and 5, and thus for the position of the cuts, as functions of wi, W2 and ai\j. In 
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particular we get: 

1 1 cos -KWi — COS TXW2 



X 



7r^(Q:2|2/2)(cOS^ T{W2 + COSTIW2 COSTTZWi — 2) + 7ra2|l S\miW2{cOS -KWi — COS7rZi;2) 



/i*-^-* = , sin 7ri(;2 (cos TTiyi — cos7r?x;2) 

V4 - n2 1 - 

^7r^a2|2/2 

X 



7r2(a2|2/2)(cos2 7rzi;2 + cos 1^)2 cos ttzwi — 2) + T^a2\\ sin 7ri(;2 (cos ttici — cos txw2) 



d) a/4 - -7r^(ai|2/2)cotan(7rtt;i) - Tram 
7+ = ■ 
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2 A/7r^ai|2/2 

(1) ^ 



v^4 — 7r^(Q;i|2/2)(2 — costtzwi) + ttciim simiwi 



2 sin TTWi ■\/7r^ai|2/2 

(5.24) 

Since si(a:(f )) is also a rational function of cos tiv with a pole at f = W2, it is possible to 
derive from Eqns. ( 4.12|) -( ^.13| ) expressions for h^^\h^'^\ together with the coefficients 
(A, 5) defined through 

s^{x{v)) = \-~5 '^"'"'''"-^ , (5.25) 

cos TXV — COS 'KW2 

in terms of Wi, W2 and ai\j. In this way, we find the alternative expressions 

h^^^ = —^=^= ^ (cos Txwi — COS 'KW2) sin nwi 

yjA - (1 - a^) 

^ ^7r^ai|2/2 

— 7r2(Q;i|2/2)(cOS^ TVWi + COS TVWi cc 



— 7r2(Q;i|2/2)(cos^ TTiyi + cosnwi cosnw2 — 2) + nai\i sin7rt(;i(cos7rii;i — cosnw2) 

1 1 cos Tvwi — cos 7rti;2 

2a/4 — 77.2 a(l — a^) sin^ ttzwi sin7rw2 
— 7r^(Q!i|2/2)(cos^ TTWi + COS TTWi COS 7rw2 " 2) + Train sin7rt«i(cos7ri(;i — cos7ri(;2) 



7r2(ai|2/2)^7r^a2|2/2 



r2) v/4 — 7r (Q;2|2/2)cotan(7rt(;2) + 7rQ;2|i 
7+ = ■ 
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2 A/7r^a2|2/2 

(2) ^ 



2) a/4 - 7r^(a2|2/2)(2 + cos7rzx;2) - 7ra2|i sin7rw2 



2 sin7rzi;2A/7r2(a2|2/2) 

(5.26) 

Let us stress that imposing Eqn. (|5.7| ) is then fully equivalent to demanding that the 
two above expressions (|5.24|) and (|5.26|) for h^^^ (or for h^"^^) are indeed equal. 
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Fig. 9: The function W2 = 4>b;r{wi) for h = 0.3 and several values of 
r. For r = 1, we have 4'b;r=i{wi) = 1 — wi (red line). Above this line, 
we have plotted 4)b;r for r — 3/2, 2, 3, 5, 7, 10 and 20. Under this line, we 
find (pb-r for the inverse values (2/3,1/2, etc.) by axial symmetry. The 
lower right quadrant corresponds to Wi €]0, 1[ and W2 e]0, 1[. In the lower 
left quadrant, wi becomes pure imaginary, more precisely in iR!j_, while 
W2 remains in ]0, 1[. In the upper right quadrant, Wi remains in ]0, 1[ 
while (1 — W2) belongs to iRi^. Finally, in the upper right quadrant, wi 
and (1 — W2) both belong to iRi^. Notice that (f)b;r is continuous across 
the axes wi = and (1 — W2) = 0. Once r is fixed, hence {wi,W2) lies 
on one of the lines in the graphic, the value of a selects one particular 
point on this line according to ( ^.7|) . For instance, a = 1 corresponds to 
points at the intersection of the fixed r-lines and the dashed magenta line 
wi = W2, increasing a corresponds to intersection points moving towards 
the lower right corner. Furthermore, the positivity conditions exclude a 
domain (here plotted qualitatively in green) in this lower right quadrant, 
without affecting other the quadrants. This leads to the existence of some 
maximal value (tyi)max[?'] admissible on the non-generic critical line, which 
in turn fixes a maximal value aniax[?']- Changing the value of b does not 
modify qualitatively (and even quantitatively when b is between 0.1 and 
0.4) the graphs above, except for 6 —t- 1/2 (i.e. n — >■ 0) where the forbidden 
green domain is pushed towards the axes. In other words, the positivity 
condition is satisfied everywhere at b = 1/2. 
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Now, we can use the set of functions introduced in Eqn. ( 4.47 ) together with the 
differential operators of Eqn. ([4.49| ), and obtain the Laurent coefficients a^j by solving 
the 4x4 system ensuring that the densities vanish at the edges. The system itself reads 



Vi [ cos rrbwij + V2 [ cos Tvb{l — W2)] = 

r'i[cos7r(l - - P2 [cos7r(l - - zws)] =0 

"sin7r6(l — 



sin7r(l — b)wi 
sin nwi 



"sin7r6wi" 


-V2 




- sinTTzwi . 





+ V2 



sin7r(l — W2) 
sin7r(l — — W2] 
sin7r(l — W2) 



0, 
0, 



(5.27) 



and its solution yields cumbersome expressions for ai\j in terms of wi, W2, u, r. We may 
plug these expressions back in the extra relation, which can be written here as 



Vi[cosn{l + b)wi] - V2[cosn{l + b){l - W2)] =0. 
The result is that Eqn. ( |5.28| ) becomes equivalent to a relation of the form 

Kbiwi, 1 - 1^2) 

r = 

Kb(l - W2,Wi)' 



(5.28) 



(5.29) 



where Kb{w,w') is a complicated trigonometric function which we will not copy here. 
We find graphically that, for a given r > and wi g]0, l[UiR;^, only two solutions 
for W2 exist in ]0, 1[U (1 — iR!^). More precisely, there exists two continuous functions 



■'b:r. 



k=i,2 such that (tfi, (f>l^}{wi)) is a solution of Eqn. (|5.29| ). However, we find that 

only one of them can yield a spectral density p*-^^ which is locally positive near 7^^''. 
In this case, local positivity implies that wi < (tyi)maxM or wi G iR!^. We call this 
function (pb-r- Numerically, it does not depend very much on b far from or 1/2. We 
give its plot for b = 0.3 and several values of r in Fig. 9. If r = 1, (j)b;i{wi) = 1 — wi 
is an obvious solution, which indeed satisfies the local positivity constraint, provided 
wi < (tyi)max[?' = 1] or wi G IR!^. We shall return to this case in Section 5.5 below. 

To summarize, Eqn. ( ^.29|) and Eqn. ( ^.7| ) (or equivalently the matching of the two 
expressions for h^^^ in ( |5.24|) and ( |5.26| )) provide two equations determining entirely wi 
and W2 as a function of r and a. For illustration, we have plotted in Fig. 10- (a) the 
values of h^^^ and h^"^^ at b = 0.4, u = 1 and for varying r (i.e. u^^^ = -\/r = l/w^^^) for 
the particular value a = 2. Surprinsingly, we see that in this case h^"^^ > h^^^ whenever 
r > 1, in contrast with what we had for a = 1. There exists therefore an intermediate 
value of a where this "exchange" takes place. For b = 0.4, we find that this exchange 
occurs for a ~ 1.7, as displayed in Fig. 10- (b). 

• The results simplify drastically in the limit n — )■ 0. We find in this case the two 
conditions 



sin^(7r 1^2/2) 
cos"^(7r 101/2) 



4a^(l-|-cos7rwi)(l — cos7rw2) =4 — (costtwi -|-cos7rt(;2)'^, (5.30) 



36 



a = 2 



a = l.l 





(l)/„(2) 



(a) 



(b) 



Fig. 10: Critical values of h^^^ and /i*-^-* as a function of r = u^"^^ /u^"^^ for 
several values of a, and the values h = 0.4, u = 1. In (b), we have chosen 
the value a = 1.7 around which an exchange between h^^^ and h^"^^ occurs 
(compare Fig. 7 and (a) here). 

which fix If; 1 and W2 in terms of r and a. More explicitly, we find 

2 2 r2/3 

1 + cos(7ri(;i) = - — , 1 — cos(7rw2) 



Va, 



where 



Va^r = \JAa?'r'^/^ + (1 - r2/3)2. 
This leads eventually to the following parametrization of the critical variety 



(5.31) 



(5.32) 



,,-l/2„-7/12 ,^ 1_Lr2/3xl/2 



(r7„,, - 1 - r^/^) 



4(a2-l) '■'"''^ ^ ' ^Vr7,,, + l-r2/3 



,,-l/2„-l/12 ,1 ^2/3.1/2 

^(2) _ ^ ' ^ ' 1 .„2/3\ /^ci,r- + J- — ^ \ ' 



4(a2 - 1) 
7l^^ = ±2z/i/2ri/4, 
7l') = ±2^.1/2 



?7a,r - 1 + 



(5.33) 



Here, it turns out that (tyi)max[?'] = 1, i-e. one can check explicitly that the spectral 
densities are locally positive at their edges for all positive values of a, r (in there words, 
the forbidden green domain in Fig. 9 is pushed toward the axes when n — t- 0). For 
illustration, we have plotted in Fig. 11 the values of h^^^ and /i*-^-' at n = 0, w = 1 and 
for varying r (i.e. u^^^ = i/r = l/w*^^)) for the particular values a = 1, a = 2 and 
a = 4/3 at which an exchange between h'^^^ and /i*-^-' occurs. Notice that, in the limit 
a — )■ 1, we recover the results of Eqn. ( [5. 20] ). Notice also that 7*^*-* = ±2a/ is totally 
expected as n — )■ corresponds to having a single loop separating a red and a green 
domain, each reduced to a tree, with respective weight per vertex u^^^ and m'^^^. 
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Fig. 11: Critical values of h^^^ and h^'^^ as a function of r = u^^^ /u^"^^ 
for several values of a, in the limit n — )■ 0, and for u = 1. The exchange 
between h^^^ and h^"^^ when passing from (a) to (b) occurs at a = 4/3 
(plotted in (c)). 

• We point out that, at n = strictly, we have 6 = 6=1/2, and therefore the dominant 
contribution to the spectral density comes from the sum of the coefficients in front of 
(j^{v ± Wj) and (h{v ±Wj). So, we need this sum to be positive, rather than cancellation 
of the coefficient of i '"'j) alone. In other words, the extra relation introduced 
in Section 5.2 is no longer necessary, and the critical variety has codimension 1. Its 
equation is obtained by demanding only si(7^^'') = 7^^'', namely 

1 - 2a^A^ h^^^ - 2a^A^ h^^'^ + 4{a^ - 1) Vm(i)w(2) = 0, (5.34) 

and the positivity condition is satisfied for all positive values of the parameters. 
Eqn. ( ^.33|) is simply the point on this line which satisfies the (no longer required) 
extra relation. From a combinatorial point of view, we do know if this point has a 
special meaning. 

5. 5. The symmetric fully-packed model 

The symmetric fully-packed model corresponds to having u'^^^ = u^"^^ = u, i.e. 
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r = 1, and h^^^ = h^"^^ = h. In this case, the condition ( 5.29 ) yields W2 = ^ — ifi, 
and the problem reduces to the 0{n) model with bending energy studied in 0, with 
moreover (7 = in the notations of that paper. We may therefore refer to the simpler 
analysis of to deduce the value of wi (hence of W2)- From the expressions of 0, we 
find that it is fixed by the condition 



6sin(27rwi)(2 + cos 2Txhwi) — sin(27r6?i;i)(2 + cos 2'kwi) 
siv? Tiwi {b sin 2nwi — sin 27tbwi) 



(5.35) 



Note that wi G [0, 1] for a > |(1 + 6^) (with in particular wi = 1/2 for a — 1), while 
Wi is purely imaginary otherwise. 



h 



. 5 



. 3 



b = OA 




a 



Fig. 12: Critical line in the [h, a) plane for the symmetric fully-packed 
model, at 6 = 0.4 and u = 1. In the blue part, wi is pure imaginary, while 
in the green part, it belongs to ]0, 1[. 



The non-generic critical point is then found at the symmetric value 



H 



\/2 — n\a'^ — l\^/u 



(5.36) 



with 



H 



Ah cos^ TTWi (—6^ sin^ ttwi cos ttwi — h sin ttwi sin 2Tcbwi + 2 cos ttwi sin^ nbwi) 

sin^ 7iwi{bs'm2Trwi — sin27r6wi) 



(5.37) 

For illustration, we have plotted the value of h versus a for b = 0.4 in Fig. 12. 

The condition that the spectral density p*-^-* be non-negative for x — )■ (7^^'')" reads 

wi < ('ii'i)max = (w^i)max[l] where (lOi)max > 1/2 is fixed by the explicit condition 
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a 



max 



lOr 



0.25 0.5 0.75 



1.25 1.5 1.75 



n 



Fig. 13: Maximal admissible value of a for non-generic critical points in 
the symmetric fully-packed model, as a function of n. Cmax ranges between 
+00 (at n = 0) to 2 (at n = 2), and assumes the value 4 at n = 1. 

/b[(Wl)max] = 0, with 

fb{w) = cos7r(6 — b)w — cos7r(4 — b)w + 2 cos7r(2 — b)w + 2 cosrcbw 

— cos 7r(2 + b)w] + 6[cos 7r(6 — b)w + cos 7r(4 — b)w — 4 cos 7r(2 — b)w — 4 cos irw 
+ 3 cos 7r(2 + b)w + cos 7r(6 + 4)zi;] +8 cos ttw sin irbw sin 7r(3 — 2b)w. 

(5.38) 

From ( ^.35|) , this condition amounts to the condition a < Cmax with Omax given by 
( p.35| ) with wi = {wi)uiax- We have plotted in Fig. 13 the value of amax as a function of 
n. It decreases from Cmax = C)0 for n = (6 = 1/2) to a^nax = 2 for n = 2 (i.e. 6 = 0). 

Note that, technically, when a < a^ax: the dominant order of the spectral densities 
when V — )■ — ioo is given by the terms involving (i,{v ± Wj). Since we obtain locally 
negative densities when a > Omax, the coefficient of this (generically) dominant order 
should vanish at a = Omax, yielding new exponents different from those of the dense 
phase. For the symmetric fully-packed model at a = amax, we expect that the new 
dominant order is given by the first subleading terms involving (h{v ± Wj), recovering 
in this way exponents of the dilute phase of the 0{n) model. 

• Forn = (6 = 1/2), Eqns. (|05|)-(^37D reduce to 



1 



a 



H 



cos nwi 



l-fCOSTTlWi 2 \1 -f- COSTriWi 

in agreement with (|5.31[ )-( |5.33| ) for r = 1. We have then [wi) 
• Forn = 1 (6 = 1/3), Eqns. (|05|)-(p:^ reduce to 



(5.39) 



1 and ttr 



oo. 



a = 



(l + 2cos(^ 



0) 



2 ' 



H 



2 (l-2cos(^))^ (5 + 2cos(^)) 
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(l + 2cos(^))^ 



(5.40) 
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We have then (lyi)max = 3/4 and Omax = 4. Let us now show that the disappearance 
of a non-generic critical solution beyond a = amax may be associated to a spontaneous 
domain symmetry breaking in the model. As already mentioned in Section 2.2, the 
n — 1 fully-packed loop model with curvature weight a is equivalent to an Ising model 
on random tetravalent maps with an Ising coupling of the form 1 + (a — l)5<j<j' and a 
magnetic-like coupling of the form + (/i*^^'')^5<j,_i which, for h^^^ = h^"^^ = h 

reduces to a weight per vertex of the tetravalent map. Up to a global factor w^, we 
may also transfer the weight u per green or red vertex of the original triangulation, i.e. 
per face of the tetravalent map, to the vertices of the tetravalent map (since a planar 
tetravalent map has two more faces than vertices), resulting in a total weight h'^u per 
vertex of the tetravalent map. 

Now we may use the known solution for the above Ising model [^5[ (see also ||26|| 



and setting H = (1 — a'^)^/i^M (in agreement with ( ^.36|) at n = 1), the solution 

of the Ising model may be characterized by the equation (see for instance Eqn. (6.4) of 
with the correspondence g ^ h'^u, gz ^ a, x H) 

P = l + 3H'P' + a'^Y^^pP' (^-^l) 

where (1 — a^)P is some appropriate generating function of the model (whose precise 
interpretation is irrelevant to the discussion here). The locus of critical points is obtained 
by differentiating this equation with respect to P, hence by demanding that 

After factorization, this yields two branches of solutions: a — {1 — 3HP)^ and (1 + 
3HP) = 0, corresponding respectively to a non- generic and a generic critical point. 
Plugging these values back in ( p.41|) and eliminating P yields the following critical 
values 

2 

H = -(1 — 3a -f 2a^/^) non-generic critical. 
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2 « (5.43) 



3a „ .... 

H = generic critical. 

36 ^ 

The passage from the non-generic to the generic critical points occurs at the value 
a = 4 (with a generic critical point when a > 4) where the two determinations of H 
are identical (note that these two determinations match up to order (a — 4)^ in the 
vicinity of a = 4). In terms of spin variables, this passage corresponds precisely to an 
Ising transition, with a spontaneous breaking of symmetry between the and —1 spin 
domains for a > 4. 

We see that the parametrization ( ^.40|) reproduces precisely the value of the non- 
generic critical point while the value amax = 4 corresponds precisely to the value of 
the Ising coupling at the Ising transition. We may thus associate the disappearance of 
the non-generic solution for a > a^ax with a spontaneous symmetry breaking between 
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red and green domains in the initially symmetric twofold fully-packed loop model. We 
expect that this conclusion, obtained for n = 1 via the above equivalence, also holds 
when n 7^ 1. Heuristically, for a large enough curvature weight, the loops are forced to 
encircle small domain of one color in a background of the other color. Note that, for 
n 7^ 1, we have however no direct proof of the symmetry breaking. 



6. Conclusion 

In this paper, we have shown that the nested loop approach, initiated in , can be 
applied to more complicated loop models than those of the 0(n)-type, which incorporate 
a bending energy a for the loops (as in [^), and possibly display a domain symmetry 
breaking. With this method, one can obtain a simple and elegant combinatorial deriva- 
tion of the functional relations satisfied by the generating functions of the model on 
planar maps, which may in some cases be solved analytically. This framework includes 
the Q = n? Potts model on general planar maps, which is equivalent to a particular 
fully-packed loop model. 

We find that imposing an explicit domain symmetry breaking (r 7^ 1) does not 
destroy the presence of non-generic critical points: the non-generic critical variety has 
codimension 2 instead of 1. More precisely in the fully-packed case, for each value of 
(r, a) with a < amaxM, there exists a non-generic critical point (/ii,/i2)- We underline 
that, whenever r 7^ 1, h^^^ 7^ h^"^^ so that the critical point is not equivalent to that 
of a dense 0{n) model, though we find that they lie in the same universality class. 
This conclusion applies to the critical Potts model on general random maps, which is 
not self-dual whenever Q 7^ 1, and does not map to a critical 0{n) model. Even with 
fully-packed loops, we expect that one may reach the universality class of a dilute 0{n) 
model by simply tuning the curvature weight to a = Omax [r]. We expect that a^a.^[r\ is 
always strictly larger than 1 and therefore the dilute phase cannot be reached strictly 
speaking in the Potts model (which assumes a = 1) without defects. For a > ttmaxf^], we 
expect a phase separation between regions with high concentration of red/low concen- 
tration of green, and vice versa. For the symmetric case (r = 1), this phase separation 
corresponds to a spontaneous breaking of symmetry between red and green domains. 
This is corroborated by known results for n = 1, as the model is then equivalent to a 
Ising model on tetravalent maps. 

Others statistical physics model on planar random maps, such as the 6-vertex 
model solved in [^-|j2^, the ADE models and generalizations thereof, might be 



treated with the same combinatorial formalism. The possible solution of those type of 
enumeration problems in all topologies by a "topological recursion formula" is currently 
under investigation ||31|| . 
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